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Abstract 



The chiral and scale anomalies of a very general class of non local Dirac 
operators are computed using the (^-function definition of the fermionic deter- 
minant. For the axial anomaly all new terms introduced by the non locality 
are shown to be removable by counter terms and such counter terms are also 
explicitly computed. It is verified that the non local Dirac operators have the 
standard minimal anomaly in Bardeen's form. 
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I. INTRODUCTION 



Local field theories provide the traditional setup where the implementation of space-time 
symmetries becomes rather simple. On the other hand, effective theories are not necessarily 
local, although an appropriate choice of degrees of freedom can make them local [|TJ. An 
outstanding example is QCD in the domain of low hadronic energies, where light quarks and 
gluons are dressed by the interaction making the effective theory look highly nonlocal || in 
terms of these degrees of freedom and hence a sort of dynamical perturbation theory would be 
needed [|[. This has produced a wealth of work mainly based on Dyson- Schwinger equations 
properly constrained by the relevant Ward and Slavnov- Taylor identities || . In this respect 
anomalies provide an interesting playground to study the interplay between low and high 
energies both in the local and in the nonlocal case. They are triggered by the ultraviolet 
regulators which unavoidably violate some classical symmetries but their physical effect is 
formulated as a low energy theorem. The question whether or not the nonlocal interaction 
can be implemented without spoiling the anomaly has been previously discussed 0-0] for 
some specific processes like e.g. tt° — > 27, 7 —>■ 3tt and 2K — > 3tt and regarding the chiral 
anomaly. In this paper we study the question for all processes and both chiral and scale 
anomalies. Rather than computing specific processes one by one we just prove that the new 
terms generated by the non locality can be subtracted by adding suitable counterterms. 



II. NON LOCAL DIRAC OPERATORS 

We will consider Dirac fermions in the flat D-dimensional Euclidean space-time R D 
endowed with internal degrees of freedom collectively referred to as "flavor". The class of 
Dirac operators to be considered here is 

D = D L + M. (1) 

The term D^, the local component of D, is a standard Dirac operator 

D L = 7^ + Y (2) 

The Dirac gamma matrices are anti-Hermitian (we will follow the conventions of ||), P M = 
idfj, and Y is an arbitrary matrix-valued function in flavor and Dirac spaces. It will be 
convenient to regard Y as a function of the position operators X^, defined by X^i/j^x) = 
x^ip^x), so that Y is a multiplicative operator in the Hilbert space of fermions. 

(Ytfj)(x) = Y(x)tfj(x) . (3) 

The term M is a purely non local, more precisely bilocal, operator also with arbitrary 
structure in flavor and Dirac spaces, 

(MV0(*) = Jd D yM(x,y)ij(y). (4) 

By purely non local we mean that M is softer in the ultraviolet sector than any multiplicative 
operator, that is, the distribution M(x, y) is less singular than the Dirac delta 8{x—y). More 
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restrictive assumptions on M will be made below. Further restrictions on the form of Y and 
M come from imposing hermiticity of the associated Hamiltonian in Minkowski space (that 
is, the hermiticity of J°D). In even dimensions, the Euclidean Dirac operator can be split 
into the components with and without 75, D_ and D + respectively, then unitarity requires 
D± = ±75-D±75. In the odd dimensional case, D + and D- corresponds to an even or odd 
number of Dirac matrices, respectively, and the unitarity condition becomes D*± = ±D±. 

Many of the concepts used for standard local Dirac operators apply directly to the non 
local case. We define a symmetry as any transformation of ip(x) and ip(x) that can be 
compensated by a corresponding transformations of the external fields Y and M (within 
the class of fields considered) so that the action / d D xipDip remains invariant. Presently, we 
will consider chiral transformations in two and four dimensions. Scale transformations will 
be treated in section |V[ Chiral transformations are defined as in the local case, namely, 

where a(x) and j3(x) are Hermitian matrices in flavor space only, regarded as multiplicative 
operators on the fermionic wave functions. The particular cases a = and j3 = correspond 
to vector and axial transformations, respectively. In the infinitesimal case 

SD = 5 V D + 6 A D = [i/3, D] - {ia l5 , D} . (6) 

Because the chiral transformations are local, both Dl and M transform covariantly sepa- 
rately, that is, 

5D L = D L ] - {m 7 5, D L } , 6M = [0, M] - {iaj 5 , M} . (7) 

Note that the bilocal structure of M implies that local factors at each side of the operator 
are taken at different points, i.e. M(x,x') -> e i(i{x) - ia{xh5 M(x, x ') e -i^')-^')^ . The 
effective action of the fermions in presence of the external fields Y and M is defined as in 
the local case, namely 

W(D) = -log J VtpVip exp |- J d D xip(x)Dij(x)^ = -TrlogD. (8) 

Here, Tr stands for trace over all degrees of freedom and some renormalization of the ul- 
traviolet divergences is understood. The (consistent) anomaly is defined as the variation of 
the effective action under infinitesimal chiral transformations. Since we will be considering 
a (^-function renormalization of W, there will be no vector anomaly, 

5 V W = 0, 8 A W = A A ■ (9) 

Correspondingly, the same current conservation formulas valid for the local case can be 
written here, 



= J d D x$(x)[i{3,D]i>(x)) Q (10) 
-A A = J d D x (4j{x){ia-f 5 , D}4j{x)) q . (11) 
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(The symbol ( )q stands for quantum vacuum expectation value.) In particular the term 
7 M P At in D in the right-hand side yields, after integration by parts, the divergence of the 
fermionic vector and axial currents whereas the other terms in D, local and non local, 
represent the explicit chiral symmetry breaking due to the external fields. On the other 
hand, the left-hand side shows the anomalous breaking of the axial current conservation. As 
will be shown below, the effective action can be renormalized so that only the local fields 
Y contribute to the anomaly, and moreover, only the standard minimal Bardeen's axial 
anomaly needs to be retained. 

For the purpose of doing detailed calculations we will assume that the non local operator 
M admits an expansion in inverse powers of P p for large P p of the form 

M = + M pv ^f + M, v p^ + ..- (12) 

The coefficients M Pl ,,, Pn are multiplicative operators and they are completely symmetric 
under permutation of indices. For convenience the P M has been put at the right. It should 
be noted that this choice does not exhaust all possible non local operators. For instance, for 
each given fc = 0,l,2,..., the class of operators 

M « = EEM-:;f 5 ^^ 7 ^ (is) 

includes all classes with lower index k as particular cases. Our choice k = is the simplest 
one but it is still non trivial and enjoys the essential property of being closed under chiral 
transformations []. The better way to obtain the transformation properties of M is by 
introducing a family of operators associated to M as 

M(p) = e ipX Me-' ipX (14) 

where the momentum p^ is just a constant c-number. Effectively, M(p) corresponds to 
make the replacement P^ — > P p + p^ in M. The function M(p) admits an expansion in 
inverse powers of p p similar to that in eq. (|12D , namely 

M(p)=M/f + M M ^ + ---. (15) 

The two lowest coefficients are given by 

M p = M M , 

= + t^ pa M p P a , (16) 

where we have introduced t pupa = 5 pu 5 pa — 5 pp 5 va — S pa 5 up . It should be noted that the 
coefficients M pi . are not multiplicative operators. One useful property of M{p) is that 



1 This choice is also realistic since it accommodates the operator product expansion estimate of the quark 
self-energy £(p 2 ) ^pi^oc (logp 2 ) d_1 / 'p 2 with d the anomalous dimension of the quark condensate 'tpip ( see 
ref. §) 
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it transforms covariantly under chiral transformations. Indeed, if Mq = Vt\MVt2 for two 
multiplicative operators fii )2 , 

Mn(p) = e ipX M n e~' ipX = ^ x M(p)0 2 . (17) 

As a consequence, the coefficients are also chiral covariant 

SM Mn = [i/3, M W „J - {za 75 , M W ...J . (18) 

From here it is immediate to derive the transformation of the original coefficients Aff M1 .. <Aln . 
For the two lowest order coefficients one finds 

5M p = M„] - {ia 75 , M M } , (19) 
5M„ U = [i(3, M pv \ - {iaj 5 , M pv ) + t^ pa M p (d a (3 + d a a l5 ) . 

In general, the variation of each coefficients involves those of lower order. This shows that 
the class of non local operators considered carries a representation of the chiral group. 

There is another essential requirement which is also satisfied by the particular class of 
non local operators M considered, namely, it is closed under Hermitian conjugation (see 
section |TV| ) and so the coefficients can be chosen so as to satisfy the unitarity requirement 
stated above. Another remark is that the classes of operators corresponding to set to zero 
the first n coefficients in M are also closed under chiral transformations and Hermitian 
conjugation. 



III. THE AXIAL ANOMALY 

In order to compute the axial anomaly, we will adopt the ^-function renormalization 
prescription combined with an asymmetric Wigner transformation. This method, as well as 
several of its applications, is presented in great detail in || . Since the techniques required 
in the present non local case are an immediate extension of those used in that reference 
here we will emphasize only the new issues introduced by the non locality. The (^-function 
effective action is given by |To| , pT| 



W(D) = -±Tr(D*) s=0 , (20) 

where s = is to be understood as an analytical extension on s from the ultraviolet con- 
vergent region Re(s) < —D. The key point is that for sufficiently negative s there are 
no ultraviolet divergences and formal operations become justified. By construction, the 
^-function renormalized effective action is invariant under all symmetry transformations as- 
sociated to similarity transformations of D, thus in particular it is vector gauge invariant. 
On the other hand the axial anomaly takes the form 

A A = Tr (2ia l5 D s ) s=() . (21) 

The operator D s can be obtained from 

D° = - [ (22) 
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where the integration path T starts at — oo, follows the real negative axis, encircles the origin 
z = clockwise and goes back to — oo. Using the Wigner transformation technique ||, the 
anomaly can be written as (a similar expression holds for the effective action) 

A ^ = -Iw^dI ^ s tr(0|2*a 757 -i 10)| (23) 

J (2n) LJ Jr 2m D(p) — z 5=0 

Here tr stands for trace over Dirac and flavor degrees of freedom, |0) is the zero momentum 
state normalized as (x|0) = 1, thus P^O) = (0|P M = 0. Further 

D(p) = e ipX De- ipX =j> +D L + M(p) . (24) 

The integration over z should be performed first, since it defines the operator D s , then the 
integral over p which corresponds to take the trace over space-time degrees of freedom and 
finally s is to be analytically extended to s = 0. The simplest way to proceed is to introduce 
a mass term, i.e., to apply the formula to the Dirac operator D + m and then make an 
expansion in powers of + M(p), letting m — > at the end. In this way the following 
expression is derived 



f d D p r dz s ... {D L + M(p))((j>+z-m){D L + M(p))) 

A A => 7 TF7 / * tr 2zcry 5 ; wwn '— 0) 

fc Q J {2-n) D Jv2m x 1 ' (p 2 + (z - m) 2 ) N + l 1 1 



=0,m=0 

(25) 



This formula has been simplified using the cyclic property for the trace in Dirac space. 

From the expression it is clear that for sufficiently large iV the integrals become ultraviolet 
convergent. When this happens s can be set to zero directly and T no longer encloses any 
singularity thus the integral vanishes. Using this insight we can expand M(p) in inverse 
powers of p^ and keep only the divergent terms. Using eq. (|i~5l) , the anomaly can be written 
as a sum of monomials each given as a product of Dl and the various coefficients M^...^ 
raised to different powers. The canonical dimension of the monomial, g, is obtained noting 
that the canonical dimension of is 1 and that of M ^ ... Mn is n + 1 (the various factors 
of m, z and p^ do not count to compute g). Doing the usual dimensional analysis, it is 
easily established that after integration each monomial comes with a factor m 7 where 7 is 
the degree of divergence of the term and further 7 = D — g + s. As noted above, terms 
with a negative degree of divergence vanish. On the other hand, the terms with a positive 
degree of divergence also vanish after taking m = 0, since m is raised to a positive power. 
In summary, the only contributions which need to be retained are those logarithmically 
divergent which correspond to monomials with scale dimension g = D. In two dimensions, 
the only relevant terms are those of the form D 2 L + M^, whereas in four dimensions they 
are D L + D 2 L M ^ + D^M^ + M ^ + M ^ up . A further restriction comes from Euclidean 
rotational invariance on p^. As a consequence, it is immediate to see that the terms M ^ in 
two dimensions and M^up in four dimensions cancel. 

After an angular average over p^, the indicated integrals on p^ and z can be carried out 
directly with the integral 1\ given in ||. The result for the two dimensional anomaly is 

A A = (2ia l5 D 2 L ), (26) 
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and in four dimensions 



A A = (2ia ls (±Di + ^D L { lfl , D L } 2 D L 

+~ (M m {7m, D l }D l + D L { lfM , M,}D L + D L { lfl , D L }M„) 

The notation (/) stands for 

</>=(4^ tr WWI°>- (28) 

Observe that, even for non local Dirac operators, the anomaly is a local polynomial of 
dimension D constructed with and the external fields Y and M^,,,^. This is a general 
property of all anomalies since only ultraviolet divergent terms can contribute to them. 
This puts a restriction to the number of non local coefficients that can appear in the axial 
anomaly, namely, only coefficients of dimension at most D can be relevant. In terms of the 
kernel M(x,y) in eq. @, it implies that non local components in the Dirac operator with 
kernels which are piecewise continuous with jump discontinuities do not contribute to the 
anomaly. The detailed calculation shows that actually the coefficients with dimension D 
(M^ in two dimensions and M in four dimensions) do no have a contribution either. In 
particular, in two dimensions there is no non local contribution to the axial anomaly. 

The expressions found for the anomaly can be put in a more usual form, in terms of 
vector and axial fields, scalar fields, etc, by making two observations. First, after taking 
the Dirac trace, the operators that appear there are actually multiplicative, that is, all P M 
appear inside commutators only. The simplest way to see this is by formally replacing every 
P M by P M + a M where a M is a constant c-number, and checking that all a M -dependence cancels. 
Second, for a multiplicative operator f(X), 

</(*)> = J dPxtrf{x) . (29) 

Since the regularization preserves vector gauge invariance, the axial anomaly is also 
invariant. In our expression for the anomaly, this is a direct consequence of the operators 
there being multiplicative. Indeed, any operator / constructed with the gauge covariant 
blocks .Dl and M^,,,^ is also covariant, i.e., / — > fl/Q -1 . If in addition / is multiplicative 
(/) is invariant. Note that ( ) is not a trace and so the cyclic property does not hold for 
arbitrary non multiplicative operators. 



IV. ESSENTIAL AXIAL ANOMALY 



Presumably due to its topological connection ]nj, the axial anomaly is a very robust 
quantity. It is not affected by higher order radiative corrections , and remains unchanged 
at finite temperature and density . It gets no contributions from scalar and pseudo scalar 
fields |15| , tensor fields |T6|Jl7| J^] or internal gauge fields, i.e, transforming homogeneously 
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under gauge transformations |18|,|l9|. In all known cases, the anomaly only affects the 



imaginary part of the effective action and only involves vector and axial fields. The counter 
terms can always be chosen so that the axial anomaly adopts the minimal or Bardeen's 
form [15[]. Not surprisingly, the new terms introduced in the anomaly by the non local 
component of the Dirac operator are also unessential, that is, they can be removed by 
adding a suitable local and polynomial counter term to the effective action. In other words, 
all new terms can be derived as the axial variation of an action which is a polynomial 
constructed with the external fields Y and M 



Ml •••Mr, 



and their derivatives. The canonical 
dimension of the polynomial can be at most D. 

The general proof that the anomaly can always be brought to Bardeen's form is as 
follows. Let Yi(x) denote the various external fields which specify the Dirac operator. 
Under a variation of them 



5D = 5Y,, 



3D 

dY~. 



(30) 



For convenience, we consider Y \ as multiplicative operators and put them at the right. The 
operators dD/dYi are simply constants, in the sense that they do not depend on Yj, and 
are just numbers or matrices if Y i refers to a local degree of freedom and contain when 
Y i is related to a non local term of D. Each external field defines a consistent current 
through the variation of the effective action 



5W = J d D xJ2^ (6Yi(x)Ji(x)) 



(31) 



At a formal level, the currents would be just Ji(x) = —(x\(dD/dY i )D~ 1 \x). However 
this expression is ultraviolet divergent. Within the (^-function prescription the renormalized 
consistent currents are given by 



It is also possible to define the currents in a chirally covariant manner 
definition PJJ corresponds to take the finite part, as e goes to zero, of 
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[x 



dX{x\^rD^ DDf \x) 

oYi 



(32) 
One such 

(33) 



(Recall that a(x) is Hermitian and so transform as D~ x under axial transformations.) 
Due to the presence of an essential axial anomaly, the consistent and covariant currents can- 
not (all of them) coincide. If the covariant currents were consistent they could be integrated 
to yield a chiral invariant effective action [21||. However, because both set of currents cor- 



respond to same formal definition, they must coincide in their ultraviolet convergent terms, 
and thus they can only differ by a local polynomial of dimension at most D — I 



Ji(x) = J c i(x) + Q^x) 



(34) 



Here, J \ is consistent, 3\ is covariant and Q i is a polynomial. Note that the arguments used 
to reach this relation hold in particular for the class of non local Dirac operators considered 
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in this work. This relation is already sufficient to show that the essential anomaly only 
contains vector and axial fields [§]. Indeed, let us separate the Dirac operator D into two 
components Dq and N both transforming covariantly under chiral transformations and such 
that Dq contains the term 7^-P^- Consequently, the vector and axial fields should also be 
in Dq since they mix with 7^-P^ under chiral transformations. Also note that Dq is a valid 
Dirac operator whereas N by itself is not. The change in the effective action due to passing 
from Dq to D = Dq + N can be obtained by integrating the consistent current along 
the path D t = Dq + tN with < t < 1. The result does not depend on the particular 
interpolating path since integrability conditions are satisfied. The contribution from the 
covariant current will be invariant and does not change the anomaly. The contribution from 
the polynomial current will be a polynomial. This latter term is responsible for the change 
in the anomaly introduced by adding N. Since this change derives from a polynomial it is 
removable by counter terms. Therefore taking Dq = 7 M (P A , + V fl + yields an anomaly 

involving only vector and axial fields. This anomaly will be minimal by using a chiral 
covariant renormalization for the real part of the effective action, such as the (^-function 
prescription applied to — |tr log(DD'). 

The actual construction of the counter terms can be done using the method in ref. ||. 
In order to keep the reasoning straight we will skip some technicalities and use a rather 
symbolic notation. The current can be written as J = dW/dD. Under vector and axial 
transformations, the covariant current transforms as J c — > e - * /3 «J c e* /3 and J c — > e*" 75 J c e Lai5 
respectively so that (NJ C ) remains invariant in both cases. (Recall that N transforms as 
e l ^Ne~ t/3 and e~ tai5 Ne~ ia ~ t5 respectively.) Infinitesimally this implies 

= 5 V J c := S V J C - J c ] , (35) 
= 6 A J c := 5 A J C - {ia 75 , J c } , (36) 

where we have introduced the covariant vector and axial variations 5y and 5 a respectively. 
On the other hand, since the polynomial current Q accounts for changes in the anomaly 
induced by changes in D, it should satisfy the following two equations 

- 8 A a 

5 V Q = 0, 5 A Q = -^. (37) 

Once these equations are solved, the counter terms needed to remove the extra contribution 
coming from N are 



W, 



ct 



(4vr) D / 2 f dt(NQ t ) , (38) 
Jo 



where Q t is the polynomial current corresponding to D t = Dq + tN. (Recall that we are 
using a schematic notation. In an actual calculation one has to distinguish each of the 
external fields which define D and their associated currents as discussed below.) 

The right hand side of the second eq. ( |37D can be computed from the known anomaly, 
eqs. (^,^). For instance, in two dimensions 

^ L {2iaih ,n L} . (39) 
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(The cyclic property has been used since it turns out to be justified in this case.) In eq. 
Q is the unknown. If the anomaly were the variation of a polynomial action, an immediate 
solution would be given by the corresponding polynomial current. Because the anomaly 
contains an essential part, such polynomial action does not exist. Remarkably, at a formal 
level there is a polynomial action from which the axial anomaly derives ||, namely, 

W = -{ l -D\) (40) 

in two dimensions and 

Wo = -<^£>1 + ^D 2 Ll ,D 2 Llll + ±Dh,D L ^ (41) 

+i( 7M M M £>i + ^D 2 L M, + 1 ,D L M,D L + m\ + {D L , M w })> (42) 

in four dimensions. In this context "formal level" means that the correct anomaly is obtained 
from Wq if the cyclic property is used. Actually, the operators involved in Wq are not 
multiplicative and so the cyclic property does not hold (it does hold in Dirac and flavor 
spaces but not for differential operators in coordinate space). Even so, the action Wq has 
a unique well-defined current Q = dWo/dD, which is obtained from (using the cyclic 
property to put all 5D together) 

6W = (4n) D / 2 (5DQ ). (43) 

Q Q is well-defined in the sense that it is unchanged under cyclic permutations of the operators 
in Wo- For instance, in two dimensions, 5Wq = —(5Dl Dl) and so 

Qo = ~^D L . (44) 

Because W formally gives the anomaly, Q Q is a solution of eq. ([37]). For instance, in two 
dimensions, SaQo is easily computed and gives precisely the right hand side of eq. fl39[). 
This solution is nevertheless formal because Q is not a multiplicative operator in general 
and Q must be multiplicative. This must be solved by subtracting to Q Q a polynomial 
current Q 1 which transforms covariantly (i.e., SaQi = 0) and such that Q = Q — Q l is 
multiplicative. Again, using the two dimensional case as an example, one can consider a 
new local Dirac operator Dl = The minus sign ensures that it is an admissible 

Dirac operator, i.e., there is a corresponding Y. The adjoint implies that it transforms as 
a covariant current under axial transformations. Then it can immediately be checked that 
the covariant polynomial current 

Qx = ~D L , (45) 

47T 



has the same non multiplicative part as Q in eq. (0). (In this case this is trivial to see. 
As noted above, the best way to make this check in general is to replace by P M + a M and 
see that Q — Q x is a^-independent.) 

To consider the non local case in detail, we introduce the currents as 
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5W = (A7r) D / 2 (5YJ L + 5M M J M + -SM^J^ + •••). (46) 

In four dimensions, the calculation proceeds by computing the non multiplicative polynomial 
formal currents 

6W = (^) 2 (5YQ° L + 5M,Ql + ^SM^Q%) . (47) 

The polynomials Q are well defined and have the correct axial transformation but they are 
non multiplicative. In order to construct the corresponding Q l , let us introduce the Dirac 
operator 

D = -D ] . (48) 

Such operator belongs to the same class as D, that is, it can be put as Dl + M where M 
has the form given in eq. (0). The corresponding coefficients are easily computed by noting 
that 

M(p) = e ipX Me~ ipX = -(M(p)) f . (49) 
And so, M^...^ = — „. For the two lowest coefficients one finds 

M, = -Ml , 

= -M\ v - t,„ pa [P p , Ml] . (50) 

The form of D is such that axially it transforms as a covariant current. The currents Q l 
are defined by suitably replacing some of the D in Q by D, so that Q 1 transforms covari- 
antly. For instance, terms of the form M^, ^^M^Dl, P p Dl or n f IJ ,D\ in Q correspond, 
respectively, to M^, •y^M^ l D L , P^D L and 7 M Z? L .D L in Q 1 . It is an exercise to check that 
Q — Qo — Qi is multiplicative. 

There is a technical subtlety in checking that Q x is actually an axial covariant current. 
This is because the coefficients iW A11 ... Atn are multiplicative but not axially covariant. There- 
fore, the covariance of a current does not directly correspond to the equation 5aJ"° ... Mn = 0. 
In order to derive the correct equations, the simplest method is to introduce the currents J 
corresponding to the covariant quantities M, 



5W = {A>k) d I 2 (5YJ l + 6M P J^ + -SM^J^ + ■■■). (51) 



They can be related to the currents J by comparing with eq. (£S^). For the two lowest orders 

J fj, J n "^/ivpa Pi/ J pa j 
J \iv = J pv + ' ' " • (52) 

Note that the currents J are covariant (up to anomaly) but not multiplicative due to the 
presence of P^. The covariant part satisfies 5 a J = 0. For the polynomial part in four 
dimensions, the relations become 
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Qfi ~ Qfi tJ'uvptjPuQ pa i 

Q,u = Q/xu • (53) 

A straightforward calculation shows that 5aQi = and so the Q Y constructed above have 
the correct axial transformation. This completes the construction of the counter terms in 
the non local case. 



V. TRACE ANOMALY 

The scale transformation ip(x) — > e~ as ( D ~ 1 ^ 2 ip(e~ as x) can be compensated by a corre- 
sponding transformation in D, namely, 

Y (x) - e- as Y(e- as x) , M^ n (x) - e~ as ^M^ n ( e - as x) . (54) 

Therefore, scale transformations are a symmetry of our class of non local Dirac operators. 
The linear operator which produces the scale transformation ip — > Qsip can be written 
as Q$ = e~ as ^ D ^ 1 ^ 2+x9 \ and thus the Dirac operator transforms as D — > Q^DQ^ 1 . 
Infinitesimally it implies 

5 S D = -a s {D - i[X^ D]) . (55) 
The corresponding trace anomaly, within the ^-function method is || 

A s = 5 s W = a s Tr(D s ) s=0 . (56) 
The calculation is entirely similar to that of the axial anomaly. In two dimensions one finds 

As = a s (D 2 L + ^{ lfl , D L } 2 + 7m M m ) , (57) 

and in four dimensions 
As = as^-Dt + ^{Dl{^ D L } 2 + { 7 „ D L }Dl{^, D L } + { 7 „ D L } 2 D 2 L ) 

+^({7m> D l } 2 { 1u , D l } 2 + ({ 7m , D l }{ 1u , D l }) 2 + { 7m , D l }{ 1u , D l } 2 {^, D l }) 
+ ^( llM MpD 2 L + ^D 2 L M, + ^D L M,D L ) 

-^{M^D 2 L + Dl^Mp + £> l { 7m , M^Dl + D Ll ^D L M^ + M^D^Dl) 

+^ uaj3 ^—(^ tJ ,M vjaDLjpDL + j li D L 'y a M v 'ypD L + j li D L 'y a D L jpM I ,) 
+L lfi M vla M p + l_{ la M^ ui lpDh } + _L 7M M ra/3 ) ) . (58) 

Where = 5^8^ + 5^5^ + 5^p5 av . The result is again a local polynomial of dimension 
D in the external fields and their derivatives. Unlike the axial case, the coefficients M M in 
two dimensions and M p Va in four dimensions do contribute to the scale anomaly. 
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Because scale and chiral transformations commute (in a properly defined sense), the 
crossed variations SsAv,a and Sy^As coincide and they vanish since the axial anomaly 
is scale invariant. Thus the scale anomaly must be chiral invariant. The vector gauge 
invariance of the previous expressions is easy to check noting that the operators inside ( ) 
are multiplicative. Axial invariance is much more involved in general. In the two dimensional 
case, it is immediate to see that ( / ~f tl M fl ) is invariant. In four dimensions it is relatively easy 
to check that the trace anomaly is axially invariant in the particular case of M ^ = 0, 
which, as noted previously, defines a class of operators invariant under chiral and scale 
t r ansf or mat ions . 

The scale anomaly is already minimal. It can be modified by adding polynomial counter 
terms of dimension smaller than D but this would add terms of the same type to the scale 
anomaly. 
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